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We use our numerical code, DRAGON, to study the implications of recent data on our knowledge 
of the propagation properties of cosmic ray nuclei in the Galaxy. We show that B/C (as well as 
N/O and C/O) data, including those recently taken by CREAM, and pip data, especially includ- 
ing recent R4MELA results, can consistently be fitted within a unique diffusion-reacceleration 
model. The requirement that light nuclei and p data are consistently reproduced within experi- 
mental uncertainties places significant limits on the main propagation parameters. In particular, 
we find the allowed range of the diffusion coefficient spectral index to be 0.3 < 6 < 0.6 at 95% 
confidence level and that Kraichnan type diffusion is favored with respect to Kolmogorov. While 
some amount of reacceleration is required, only a limited range of the Alfven velocity value 
(10 < V/i < 20 km s"') is allowed by a combined analysis of nuclear and antiproton data, which 
we perform here for the first time. If antiproton data are not used to constrain the propagation 
parameters, a larger set of models is allowed. In this case, we determine which combinations of 
the relevant parameters maximize and minimize the antiproton flux under the condition of still 
fitting light nuclei data at 95% C.L. These models may then be used to constrain a possible extra 
antiproton component arising from astrophysical or exotic sources (e.g. dark matter annihilation 
or decay). 
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1. Introduction 

The problems of origin and propagation of Cosmic Rays (CRs) in the Galaxy are long stand- 
ing questions and the combination of several different observations in a wide energy range is 
required to understand them at least partly. 

The most realistic description of CR propagation is given by diffusion models. Two main 
approaches have been developed so far: analytical (or semi-analytical) diffusion models (see 
e.g. (Ut] and ref.s therein), which solve the CR transport equation by assuming simplified distri- 
butions for the sources and the interstellar gas, and fully numerical diffusion models. Well known 
realizations of these two approaches are respectively the two-zone model (see e.g. |2,^^) and 
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DRAGON (Diffusion of cosmic RAys in Galaxy modelizatiON) f^. All these models involve in 
general a large number of parameters which need to be fixed using several types of experimental 
data. Their knowledge is crucial not only for CR physics but also for constraining or determining 
the properties of an exotic galactic component from indirect measurements. 

However, in spite of the strong efforts made on both observational and theoretical sides, most 
of these parameters are still poorly known. One of the reasons lies in the fact that best quality data 
on CR spectra (e.g. the ratios of secondary to primary nuclear species) were available mainly at 
low energy (£ < 10 GeV/n), where several competing physical processes (e.g. solar modulation, 
convection, reacceleration) are expected to affect significantly the CR spectra by an a priori 
undetermined relative amount. Furthermore, the uncertainties on the spallation cross sections 
and their effects on the propagated CR composition are still sizable at such low energies. 

On the other hand, the interpretation of high energy (£ > 10 GeV/n) CR data is, in principle, 
easier since in this range only spatial diffusion and spallation losses (the latter becoming less 
and less relevant with increasing energy) are expected to shape the CR spectra. Furthermore, 
other uncertainties related to the physics of solar modulation and to poorly known nuclear cross 
sections are reduced by considering only data at energies larger than several GeV/n. Hence, the 
study of high energy CR spectra allows in principle to constrain the plain diffusion properties of 
CR in the Galaxy, in particular the strength Do of the diffusion coefficient at a reference rigidity 
and its energy slope 6, and offers a lever arm to better understand low energy effects (see flu] for 
an interesting discussion about this issue). This possibility has been precluded for long time by 
the scarcity of observational data. 

The experimental situation however improved recently when the CREAM balloon experi- 
ment measured the spectrum of light CR nuclei and especially the boron to carbon ratio (B/C) 
up to ~ 1 TeV/n jl ill . Besides CR nuclear measurements, valuable complementary data were 
recently provided by the PAMELA satellite experiment which measured the antiproton to pro- 
ton ratio up to ~ 100 GeV with unprecedented accuracy 1 12]. Other valuable experimental data 
are expected to come from AMS-02 11 1311 which will soon be installed on board of the Interna- 
tional Space Station. As for other secondary nuclear species, antiprotons are produced by the 
spallation of primary CRs (mainly protons and Helium nuclei) in the standard scenario. There- 
fore, their spectrum may provide an independent and complementary check of the validity of CR 
propagation models and a valuable probe of an extra component which may arise, for example, 
from secondary production in the CR astrophysical sources 1,14 ,15.] and/or from dark matter 
annihilation or decay (see e.g. lfl6l. 17 . TsIITqIi V 

Whether the measured secondary/primary nuclear ratios and antiproton spectra are fully com- 
patible within the framework of a standard CR transport model is still not completely clear. In- 
deed, while a discrepancy between the parameters allowing to reproduce the B/C and the p/p 
was claimed in f?^ (see also fs*]), a good concordance was found in other analyses |fl6i l2l[l . 
Furthermore, even the interpretation of nuclear data alone is still confused: analyses based on 
the leaky-box and semi-analytical diffusion models favor values of 6 significantly larger than 
the ones found with the numerical GALPROP package. The comparison of such results is not 
straightforward due to a number of different assumptions. Hence, an independent analysis ac- 
counting for most recent available data is timely. 

In this work we use DRAGON |9] to constrain the main diffusion parameters against updated 
experimental data in the energy range 1 < £ < 10^ GeV/n. This code reproduces the results of 
the well known GALPROP under the same conditions. Furthermore, it allows to test the effects of 
a spatially varying diffusion coefficient. Here we use the optimized and updated version of this 
code, which now accounts for ionization and Coulomb energy losses, diffusive reacceleration 
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and convection, and exploits the performances of modern computer clusters to scan a rather 
large range of parameters under realistic physical conditions. These upgrades allow to constrain 
the diffusion coefficient normalization and spectral index, as well as the Alfven velocity va, 
with unprecedented accuracy by means of a statistical analysis of the agreement between model 
predictions and CR data including recent nuclear and antiproton data. 

In the following we will present the results of this analysis. In Sec.|2]we briefly review the 
framework of CR propagation we adopt. In Sec. [3] we describe our analysis and constrain the 
diffusion parameters. In the same section we also discuss how much the secondary antiproton 
spectrum is allowed to vary under the request that the predicted B/C, N/O and C/O ratios are com- 
patible with the experimental data. In Sec|4]we introduce an effective diffusion-reacceleration 
model which allows to match all relevant experimental data down to £ ~ 0.1 GeV/n. Finally 
in Sec.s|5]and|6]we compare them with results from other groups and discuss differences and 
implications for exotic source searches. Section |6] is further devoted to our final remarks and 
conclusions. 



2. The CR propagation framework 

Galactic CRs propagate diffusively in the irregular component of the Galactic magnetic field 
undergoing nuclear interactions with the gas present in the InterStellar Medium (ISM). Simi- 
larly to previous treatments, we assume here that CR Galactic source, magnetic field and gas 
distributions can be approximated to be cylindrically symmetric. Under these conditions, and 
in the energy range we are interested in, CR propagation of stable nuclei obeys the well known 
transport equation (Ginzburg and Syrovatskii 12211 ') 

- - v.,z,v-.,«'.i-(.-£v..)«.-±.=.„,i-^. 

= Q'ip, '■,z) + Yj Cy6nga,s('-, z)crjiN' - c/5n^,,crin(Ek)N' . (1) 

Here N'{p, r, z) is the number density of the i-th atomic species; p is its momentum; j6 its velocity 
in units of the speed of light c; cr,„ is the total inelastic cross section onto the ISM gas, whose 
density is «gas; cr/j is the production cross-section of a nuclear species j by the fragmentation of 
the i-th one; D is the spatial diffusion coefficient; Vc is the convection velocity; /3 is the particle 
speed in units of c. The last term on the l.h.s. of Eq. ([TJ describes diffusive reacceleration of CR 
in the turbulent galactic magnetic field. In agreement to the quasi-linear theory we assume the 
diffusion coefficient in momentum space Dpp to be related to the spatial diffusion coefficient by 

the relationship (see e.g. HI]) Dpp = '^^^ — — S)^^ 1^ where va is the Alfven velocity. 

Here we assume that diffusive reacceleration takes place in the entire diffusive halo. 

Although DRAGON allows to account also for CR convection, we neglect this effect in the 
present analysis showing a posteriori that it is not necessary to consistently describe all the 
available data above 1 GeV/n (see Sec.|5]l. Hence in the following we will set vv = 0. By this we 
do not mean that CR data implies that the physical value of Vc is actually vanishing but only that 
an effective description of their propagation is possible even if convection is disregarded (see the 
discussion at the end of Sec|5]l. 

DRAGON ^ solves Eq. ([TJ numerically in the stationary limit dNi/dt = by imposing the 
foUowing boundary conditions: A^(p,^max,z) = N(p,r,Zmin) - Nip,r,Zm!ix) = 0, corresponding 
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to free escape of CRs at the outer limit of the Galaxy; a symmetry condition on the axis r - 0, 
N(p, + e,z) = N(p, - e, z) (e <K 1), due to the assumed cylindrically symmetric setup; a null 
flux condition dN/dp = on the momentum boundaries, which stems for the fact that particles at 
null momentum should not lose momentum anymore. Even though the presence of reacceleration 
can effectively invalidate this condition, by producing a net flux from low to high momenta, we 
remark that this can afi'ect only the part of the spectrum close to the momentum boundary. For 
this reason, we adopt an energy grid whose extrema are well below and well above the minimal 
and maximal energy of the data set we consider. In such a way, our results are in fact independent 
of the momentum boundary conditions we impose. The spatial limits of our simulation box are 
defined by R^ax = 20 kpc and Zmax = -Zmm- We start the spallation routine from Z = 16, having 
verified that the effect of heavier nuclei on the results of the present analysis is negligible when 
compared to other uncertainties, being below the 1 % level. 

We briefly recall below the main assumptions we make for the terms appearing in Eq. ([1). 

2.1. Spatial diffusion coefficient 

The dependence of the diffusion coefficient D on the particle rigidity p and on the distance 
from the Galactic plane z is taken to be (here we assume D to be cylindrically symmetric and 
independent on the Galactocentric radius r) 

D(p,z) = DoyS-'^^j exp{|z|/z,) , (2) 

where /3 is the particle velocity in units of the speed of light c. As shown in | 9J, a vertically grow- 
ing D is physically more realistic than a uniform one and allows to get a more regular behavior 
of the CR density at the vertical boundaries of the propagation halo with respect to the case of 
uniform diff'usion. As far as the analysis discussed in this paper is concerned, however, the sub- 
stitution of such a profile with a vertically uniform D only requires a change of the normalization 
factor Dq. Generally, the value t] - I is adopted in the related literature. This parameter, how- 
ever, is not directly constrained by independent observations and other values have been recently 
considered (see e.g. Ia^ ). We neglect here a possible dependence on the radial coordinate r, 
which was considered also in We always set Zmax = 2 x z, in Eq. (|2]i to avoid border effects, 
and po = 3 GV in the following. Finally, we assume no break in the power-law dependence of D 
on rigidity, and we checked that our results do not depend on the choice of Znvdx, but only on Zi, 
which then acts as the effective vertical size of the diffusive halo. 

2.2. Cosmic ray sources 

For the source term we assume the general form 

Qi{Et,r,z)^fs(r,z)q'o(^^] , (3) 

and impose the normalization condition fs(f0,Z0) - 1. We assume fs{r,z) to trace the SNR 
distribution as modeled in [23] on the basis of pulsar and progenitor star surveys fl^. This is 



slightly different from the radial distributions adopted in H25I1 and in |2, 4] which are based on 
pulsar surveys only. Two-zone models assume a step like dependence of fs (r, z) as function of z, 
being 1 in the Galactic disk (|z| < Zd) and outside. For each value of 6 in Eq. (|2]l we fix c,- by 
requu'ing that at very high energy (£^ » 100 GeV/n) the equality ai + 6 = ji holds, as expected 
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in a plain diffusion regime. Indeed, at such high energies reacceleration and spallation processes 
are irrelevant. Here we adopt the same spectral index (y, - j, hence a, = a) for all nuclei as 
indicated by recent experimental results 



The low energy behavior of Q is quite uncertain and several different dependencies of Q on 
the velocity fi have been considered (see e.g. [2]). In the energy range explored in this work, 
however, different choices of such behavior have negligible effects. This strengthens further the 
importance of relying on high energy data to reduce systematic uncertainties. 

The injection abundances q^^ are tuned so that the propagated, and modulated, spectra of pri- 
mary species fit the observed ones. Here we choose to normalize the source spectra of Oxygen 
and heavier nuclides to reproduce the observed spectra in CRs at £ ~ 100 GeV/n. On the other 
hand. Carbon and Nitrogen (which, together with Oxygen mostly affect the B/C) injection abun- 
dances (with respect to Oxygen) are free parameters, over which we marginalize our statistical 
variables in our analysis, in a way which we will describe in Section[3] Our data basis for Oxygen 
and heavier nuclei is constituted by ACE/CRIS data L29J. For B, C and N, besides CREAM's, we 
use experimental data provided by the HEAO-3 ||30|] and CRN jsil] satellite-based experiments. 
HEAO-3 B/C data are nicely confirmed by a recent preliminary analysis of AMS-01 data ll32|] 
which, however, we do not use in this work. 

We verified a posteriori that the observed Oxygen spectrum (see below), as well as the 
subFe/Fe ratiofl are reasonably reproduced by our best-fit model. 

For the primary proton local interstellar spectrum (LIS) we adopt Jp - 1.6x10'* (£*/! GeV)"^-^^ 



'p 

,2 „ ri^\i\-\ 



(m s sr GeV) as measured by BESS during the 1998 flight II33I1 . This spectrum also provides 
an excellent fit to AMS-01 Jsill data and, as we will show below, also to preliminary PAMELA 
proton spectrum data Issll . 

What is most important here, however, is that we assume no spectral breaks in the source 
spectrum of all nuclear species. As we will discuss in Sec. |5]this point is crucial to understand 
the difference between our results and those of some previous works. 

2.5. Nuclear cross sections 

The spallation cross sections and the spallation network are based on a compilation of ex- 
perimental data (when present) and semi-empirical energy dependent interpolation formulas as 
provided e.g. in ll36l. 37, 38] (see also GALPROP, ff\ and references therein, from which data 



and some related routines have been obtained and included in DRAGON as an external library). 

For antiprotons, the main processes responsible for their production are p - pg^^, P — Heg^s, 
He - peas and He - Hcgas reactions, plus a negligible contribution from other nuclei. Similarly 
to i20l I2III we adopt the p production cross-section calculated using the parametrization given 
in Tan & Ng lf39ll . Inelastic scattering, annihilation and tertiary p (antiprotons which have been 
inelastically scattered) are treated as in f^O*]. 

In order to test the possible dependence of our results on systematical uncertainties on those 
cross sections, we performed several DRAGON runs using also a different set of nuclear cross 
sections as determined in [Cot] (see Sec|5]). 

2.4. Target gas 

The ISM gas is composed mainly by molecular, atomic and ionized hydrogen (respectively, 
H2, HI and HII). Here we adopt the same distributions as in t5„^ . We checked that other possible 
choices do not affect significantly our final results. 



To compute these ratios, of course we extended our numerical simulations up to Z = 28. 
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Following 114 ill we take the He/H numerical fraction in the ISM to be 0. 1 1 . We neglect heavier 
nuclear species. 



2.5. Solar modulation 

We describe the effect of solar modulation on CR spectra by exploiting the widely used force- 
free approximation ||42ll . prescribing that the modulated spectrum J(Ek,Z,A) of a CR species is 
given, with respect to the Local Interstellar Spectrum (LIS) Jus(Ek,Z,A), by 



where m is the nucleon mass and <I> is the so called modulation potential. This potential is known 
to change with the solar activity with a period of 11 years. It must be stressed that the potential O 
is not a model independent quantity. Rather, for each propagation model it should be obtained by 
fitting the CR spectra at low energy. The possibility of restricting our analysis to > I GeV/n 
will reduce the systematic uncertainties associated to this unknown. Above 1 GeV/n the effects 
of modulation on the secondary/primary CR ratios used in our analysis are tiny and can safely be 
accounted for by means of the simple force free approximation. 

For protons and antiprotons we use a potential which allows to match BESS98 ||33tl, AMS- 
01 Q and PAMELA fs?) proton data even well below 1 GeV/n (see Fig. [5(c)l ). Indeed all these 
experiments took their data in a period with almost the same, almost minimal, solar activity. 
Although a more complicated and realistic treatment of solar modulation, accounting for charge 
dependent effects, and the 22 year cycle change of polarity associated to solar effects, might be 
needed when dealing with p/p ratios for Ek <l GeV/n (see e.g. |43]), we decide to work in the 
framework of the force-free field approximation and show a posteriori that the data considered 
in our analysis can naturally be described in that framework. 

3. Analysis and results 

Our goal is to constrain the main propagation parameters 6, Do, Zi and va entering Eq. (|2|l. To 
this aim, we compare to experimental data our prediction for the following physical quantities: 
the B/C, N/O, C/O ratios for I < Et < 10^ GeV/n and the p/p ratio for 1 < < 10^ GeV/n. We 
will check a posteriori that also the Oxygen, proton and antiproton absolute spectra are correctly 
reproduced by our preferred models. In order to test the relevance of low energy physics on our 
constraints of the diffusion-reacceleration parameters, we perform our analysis for three different 
values of the minimal energy E„^i„. We will then motivate the choice of the most suitable value 

of £'min ■ 

As long as the propagation halo scale height is allowed to vary within the range 2 < z, < 
6 kpc (which is what we assume here). Do and z, are practically degenerate so that our results 
depend only on the ratio Do/zt- Throughout this paper we will always express this quantity in 
units of 10^*^ cm^ s"' kpc"' . We verified a posteriori that for this range of z, values, the predicted 
"'Be/'^Be ratio, which constrains the CR propagation time hence the vertical scale height of the 
propagation region when combined with secondary/primary stable nuclei data, is consistent 
with experimental data. 



7(£,,Z,A) = 



(Ell + nif' - 




(4) 
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3.1. Light nuclei ratios 

3.1.1. Method 

We already showed |@] that in order to constrain correctly the propagation parameters on the 
basis of B/C measurements it is essential to take into proper account that the main primary parent 
species of Boror0 are also affected by propagation. This holds not only for the Nitrogen (N - '"^N 
+ '^N), which gets a significant secondary contribution, but also for Carbon and Oxygen, since 
for Ek < 100 GeV/n their spectra are shaped by spallation losses in a propagation dependent 
way. Therefore, we perform our likelihood analysis in three steps: 

1. for fixed values of the propagation parameters v^, 6, and Dq/z, we vary the C/O and N/O 
source ratios to compute the (which we call;^'^ q) of the propagated, and modulated, 
C/O and N/O ratios against experimental data in the energy range I <Ek< 10^ GeV/n; 

2. for the same fixed value of v^, we finely sample the parameter space {5, Dq/z,) by using, for 
each couple of these parameters, the C/O and N/O source ratios which minimize q; for 
each of these realizations we compute the;if^ (which we call;i^g^^) for the B/C modulated 
ratio against data with E > E^in, 

3. we repeat the same analysis for several values of to probe the effect of diffusive reac- 
celeration. For each value of va we then determine the allowed ranges of 6 and Dq/zi for 
several Confidence Levels (CL). 

In Jgt] only items (i) and (ii) were performed, for = and without accounting for CREAM 
data, not yet public at that time. 

The wide energy range covered by these recent data allows us to perform our analysis using 
three different energy intervals defined by = 1,5 and 10 GeV/n respectively and by the 
same iSmax = 1 TeV/n. 

As we already stated in the above, we do not account in our analysis for light nuclei and 
antiproton data with energy below 1 GeV/n as they are affected by poorly known low energy 
physics and are not necessary to constrain the high energy behavior of the diffusion coefficient, 
which is the main goal of this work. In Sec. |4]we will show, however, that specific models which 
fit all data even below that energy can be built, adopting diffusion coefficients allowed from our 
analysis. 

3.1.2. Results 

In Tab. [U we report the best-fit model parameters, and the relative minimal ;^'g^p's, as deter- 
mined for several values of va and E„^i„. 

First of all we notice that in the highest energy range (Zsmin =10 GeV/n) the best-fit model 
values of 6 and Dq/zi are weakly dependent on the Alfven velocity. In particular, the best fit 
values of S stays in the very narrow range 0.40 0.46 varying va from to 30 km/s. This agrees 
with the common wisdom that reacceleration is almost ineffective at such high energies (see also 
Fig.S. 

The most useful results, however, are those obtained for Ej^in = 5 GeV/n since that thresh- 
old provides the best compromise between the two opposite requirements: 1) to include in the 



^ A non negligible contribution to the '"B comes from the beta decay of '"Be, which is properly accounted for in our 
analysis. 

^Every time we refer to a.x^, we mean the;;^^ divided by the number of degrees of freedom, i.e. the so called reduced 
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Do/Zj [10'' cm' s"' kpc"'] 

Figure 1: The 68%, 95% and 99% confidence level regions of DRAGON models, computed for = 5 GeV/n are 
represented in the plane (Do/zt, S). For the 68% confidence level the corresponding value of the^^ is also shown. The red 
crosses show the best-fit position. Each row corresponds to different values of the Alfven velocity: va = 10, 20, 30 km/s 
from top to bottom. Each column corresponds to different analyses: B/C (left panels), p/p (center panels) and combined 
(right panels). 
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Table 1: Best fit parameters, and the corresponding values resulting from comparing our model predictions with 
nuclear experimental data alone (B/C analysis) and with nuclear and pjp combined data (combined statistical analysis), 
as described in text. The values corresponding to E^^^ = 5 GeV/n for the combined analysis, which are used to constraint 
our models, are reported in bold. 





B/C analysis 


joint analysis 


va [km/s] 


£min [GeV/n] 


6 




x' 


(5 




x' 




1 


0.57 


0.60 


0.38 


0.47 


0.74 


3.25 





5 


0.52 


0.65 


0.33 


0.41 


0.85 


2.04 




10 


0.46 


0.76 


0.19 


0.44 


0.82 


1.57 




1 


0.52 


0.68 


0.32 


0.49 


0.71 


1.47 


10 


5 


0.49 


0.71 


0.28 


0.41 


0.85 


1.69 




10 


0.44 


0.82 


0.20 


0.44 


0.82 


0.12 




1 


0.46 


0.76 


0.33 


0.47 


0.76 


0.94 


15 


5 


0.49 


0.73 


0.26 


0.44 


0.82 


0.12 




10 


0.44 


0.84 


0.18 


0.41 


0.98 


0.16 




1 


0.41 


0.90 


0.47 


0.47 


0.79 


2.28 


20 


5 


0.44 


0.84 


0.22 


0.44 


0.84 


0.85 




10 


0.44 


0.87 


0.20 


0.44 


0.85 


0.98 




1 


0.33 


1.20 


0.40 


0.33 


1.20 


5.84 


30 


5 


0.38 


1.06 


0.20 


0.36 


1.09 


2.47 




10 


0.41 


0.98 


0.16 


0.38 


1.04 


1.61 



analysis more experimental data and 2) to work in an energy range where propagation is as less 
as possible affected by poorly known low energy physics. For example, possible charge de- 
pendent drift effects in the solar modulation (see e.g. |43, 20]) can be safely neglected in that 
energy range. Best fit parameters and confidence level contours obtained for that value of 
are showed in Tab.[T]and in Fig.[T]respectively. 

From both we notice that all considered values of are almost equally permitted by the 
B/C analysis, and that the d - DqIzi allowed region slightly moves towards low 5's and large 
D()/zi's as Va is increased from to 30 km/s. While Kraichnan diffusion is clearly favored in the 
case of low values of va, Kolmogorov becomes favored, for va > 30 km/s. The choice among 
those model, however, is difficult in the absence of an independent estimate of va- We will show 
that the antiproton/proton data break such degeneracy. 

In Fig. |2(a)| we show the effect on the B/C ratio of varying va keeping 6 and Dq/zi fixed to 
the value (0.45, 0.8) which will be motivated below. 

3.2. Antiprotons 

The statistical analysis for the pjp ratio is rather simpler than the one for B/C. Indeed, the 
secondary p production depends, besides on Dq/zi, 6 and va, only on the source abundance 
ratio He/p. This last unknown quantity can be easily fixed by looking at the measured spectrum 
of He at Earth, which is relatively well known. Therefore, we do not need to fit the source 
abundance ratio here and can directly proceed to map thexj)/p in the {Dq/zi, S) space, for several 
Va, similarly to what described in items (ii) and (iii) of the previous subsection. 

In the second column of Fig.[T]we show the statistically allowed regions in the plane (Do/zt, 6) 
for several values of va and compare them with the corresponding regions determined from the 
light nuclei analysis (first column in the same figure). The allowed CL region is significantly 
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larger than that determined from the light nuclei analysis (due to the larger experimental errors) 
and they overlap only for some values of the Alfven velocity. In fact, it is remarkable that the va 
varying behavior of those regions is almost opposite so that not all values of va are allowed by a 
combined analysis (see Sec. I3.3l l. 

3.3. Combined analysis and constraints on the propagation parameters 

A combined analysis of light secondary/primary nuclei and antiproton/proton data can be 
performed under the working hypothesis that CR antiprotons are only of secondary origin. 

We define the combined reduced as xlomb ~ 9 (^bc + A'ap/p)- The CL regions for several 
values of va are reported in the third column of Fig.^Jand the corresponding best-fit parameters 
in Tab.[T] Again, here we use only data with E > /imin = 5 GeV/n. 

As we anticipated in the previous subsection, in general the CL region allowed by the com- 
bined analysis is smaller than the B/C one. 

Indeed, while the parameter regions constrained by the B/C and p/p data nicely overlap for 
10 < Va < 20 km s"^ models outside this range do not allow a combined fit of both data sets at 
the required level of statistical significance (higher than 95%). The fact that only a limited range 
of the Alfven velocity values are allowed is consequence of the different behavior of the B/C and 
p/p ratios with va due to the different spectral shapes of these ratios. This is a new and quite 
interesting results. 

It is reassuring to notice that the results of the analysis performed for Emin = 5 and 10 GeV/n 
are practically coincident, which makes us confident that the combined analysis performed for 
Emin - 5 GeV/n probes already the purely diffusive CR regime. It is also remarkable that the 
best fit values of 6 and Dq/z, stay almost unchanged when varying va- In particular (6, Dq/zi) - 
(0.4 - 0.45,0.8) for all allowed values of va = 10 - 20 kms"' of the Alfven velocity. This 
makes us confident that the combined analysis performed for = 5 GeV/n best probes the 
diffusion-reacceleration parameters . 

Among those considered va = 15 km s"' is the Alfven velocity value which minimizes the 

of the combined analysis, hence it gives rise to the best overlap between the light nuclei and 
the p/p confidence regions. This is also visible from Fig.s |2(a)] and |2(b)] where the B/C and the 
p/p ratios computed with (6, Do/zt) - (0.45, 0.7) are plotted for several v^'s. It is also interesting 
to notice that the dependence of the p/p ratio on va is driven by that of the proton spectrum since 
the absolute p spectrum is practically unaffected by re-acceleration (see Fig |2(c)| i. We stress that 
Fig.s |2(a)|2(d)| are given here mainly for illustrative reasons since, below few GeV/n some other 
physics needs clearly to be introduced to reproduce the B/C data (see Sec. |4|. 

Since the = 15 kms"' combined analysis CL region is the largest, it also provides the most 
conservative constraints on 6 and Do/zt- They are 0.3 < <5 < 0.6 and 0.6 <D()/zt <l at 95% CL. 

It should be kept in mind that our analysis accounts only for statistical experimental errors. 
Several systematic uncertainties, however, may affect our constraints too. Among them, sys- 
tematic errors in the experimental data, uncertainties in the Galactic gas density and hydrogen 
fraction distributions and nuclear fragmentation cross sections play a major role. A detailed dis- 
cussion of the possible impact of these uncertainties on the determination of the CR propagation 
parameters is beyond the aims of this work. A thorough analysis was recently performed in 1 441] 
showing that, if low energy data are accounted for (which requires to introduce several unknown 
parameters with respect to those considered in this work) the systematic uncertainties on the 
Do, 6 and va can be comparable, or even larger, than the statistical ones. However, the former 
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Figure 2: The B/C (panel a), and the pip (panel b) ratios, as well the antiproton (panel c) and proton (panel d) ab- 
solute spectra computed with DRAGON for S = 0.45 and Do/^r = 0.8 are plotted for several values of v'a and 
compared with the respective experimental data. Dotted, short-dashed, solid, dot-dashed, long-dashed correspond to 
va = 0, 10, 15, 20, 30 km/s respectively. Here (7 = 1 which clearly does allow to match nuclear data below 1 GeV/n. For 
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Figure 3: The p absolute spectrum is shown for = 10, 20, 30 (from the left to the right panels respectively). The upper 
and lower curves correspond to the MAX and MIN models defined as in Sec. l3.4l respectivelv. 



uncertainties are significantly smaller if one considers only a subclass of models without convec- 
tion and keeping fixed other parameters which only matters a low energies, as we do in this work. 
For example, it was shown in 1.441 that for models with vv - the effect of considering different 
cross-section sets amounts to a ~ 40 % uncertainty variation of 5 which reduces to ~ 10 % if one 
considers only the most updated cross section sets. 

We verified with DRAGON that changing the GALPROP nuclear fragmentation cross sec- 
tions with those given in ifioll produces only a marginal effect on the B/C ratio. The relative 
effect of cross section uncertainties on the antiproton/proton ratio is negligible here due to the 
high statistical errors on those data. 

3.4. Maximal and minimal antiproton spectra 

The previous results clearly favor a standard interpretation of the measured antiproton spec- 
trum in terms of purely secondary production from CR nuclei. It is still possible, however, that a 
subdominant antiproton component arises from unconventional processes. In order to constrain 
such "exotic" component(s) with experimental data, one has to compare antiproton data with the 
predictions of the theoretical models validated against CR nuclei data alone. 

For this purpose we define, for each value of \'a considered in the above, a pair of MAX and 
MIN models which maximize and minimize respectively the antiproton absolute flux integrated 
in the range 1-100 GeV under the condition to be compatible with secondary/primary light 
nuclei data down to 1 GeV/n within 95% CL. 

In Fig. [3] we show the allowed ranges of the antiproton absolute spectrum for several values 
of va- Among the models considered here the absolute MAX and MIN models are those defined 
by the parameters {6, Do/zt, va) - (0.68, 0.46, 0) and (0.30, 1.2, 30) respectively. Therefore, we 
conclude that, under the hypotheses adopted in this work, p constraints on an exotic component 
should not use, as propagation models, any model whose p background prediction is lower than 
our MIN (or larger than our MAX) model, as it would be in contrast with B/C data at 95% CL. 
Hence, the most conservative constraint, under our hypotheses, arises from the request that the 
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sum of the background p predicted by the MIN model plus the exotic p component do not exceed 
the experimental data, within some CL. 

4. A comprehensive model describing all data sets down to 0.1 GeV/n 

The aim of this section is to test the consistency of our previous results with CR data below 
few GeV/n and to identify an effective model allowing to fit all available data. 

It is evident from Fig. |2(a)| that while the best fit model obtained for t] - \ provides an 
excellent fit of experimental data above few GeV/n, below that energy it overshoots the B/C 
observations. As we discussed, such a discrepancy may be attributable to a number of effects 
which, at low energies, introduce degeneracies among the relevant parameters. For this reason, 
a statistical analysis aimed to fit those low energy parameters against presently available data 
would be hardly interpretable (see e.g. |44]) and it is beyond the aims of this work. 

Here we follow a more phenomenological approach tuning only the parameter rj (see Eq. |2]i 
which sets the dependence of the diffusion coefficient on the particle velocity (a similar approach 
was followed in |44]). Interestingly we find that the choice rf ^ -0.4 allows to match light nuclei 
as well as antiproton data well below 1 GeV/n for almost the same range of S and Dq/z, values 
found for t] - \. Indeed, we checked that the rj - -0.4 and tj - \ CL regions computed for 
Emm - 5 GeV/n almost coincide (which is not the case for Emm - 1 GeV/n. In Fig.s |4(a)| - |5(c)| 
we show as our best fit model obtained for rj - -0.4, 6 - 0.5, DqIzi - 0.7, and va - 15 km/s"' 
nicely reproduces all relevant data sets. They include also the N/O and C/O ratios (with ~ 6 % 
and ~ 75 % injection ratios respectively) as well as the absolute oxigen spectrum. 

We notice that the modified dependence of the diffusion coefficient upon rigidity, which is 
the consequence of adopting a value of different from 1, can be considered as an effective 
modelization of physics taking place at low energy, including some non-linear phenomena such 
as the dissipation of magneto-hydrodynamics (MHD) waves by their resonant interaction with 
CRs iH. Since this is the same interaction responsible for CR diffusion in the ISM, such an 
effect is unavoidable at some level. Interestingly, the value of 5 used in to fit the B/C in 
the presence of MHD wave dissipation is 0.5, which is consistent with what we found here 
(differently from what we do here, however, a break in the injection index was invoked in that 
work). 

5. Discussion and comparison witli previous results 

As we mentioned above, our numerical diffusion code DRAGON reproduces the same re- 
sults of GALPROP 101 under the same physical conditions. Our analysis and main conclusions, 
however, differ significantly from those reported in several papers based on that code. 

In order to clarify the reasons of such a discrepancy, in Fig. |4(a)| and |5(a)| we compare the 
predictions of our reference diffusion-reacceleration model {6, Dq/zi, va, ??) = (0.5, 0.7, 15, -0.4), 
which for brevity we call Kmichnan model, with those obtained using the propagation parameters 
(and source distribution) of the Kolmogorov model discussed in 1^, namely {6, Dq{A G'V)/zt, va, rf) = 
(0.33,1.45,30,1)0 For the latter combination of parameters we consider two variants, repre- 
sented by the solid/dashed red lines, which differ for the presence/absence of a break at pbreak - 



''in a spatially uniform diffusion coefficient (z, = Zmax = 4 kpc) was assumed. As we already noticed, for the 
purposes of the present analysis adopting a vertically uniform rather than varying diffusion coefficient only amounts to a 
rescahng of Dq/zi- We verified that this does not affect any other result of our analysis. 
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Figure 4: The B/C (panel a), N/0 (panel b), C/O (panel c) and the oxigen absolute (panel d) spectra computed with our 
preferred Kraichnan model (blue solid line), the Kolmogorov reference model (red solid line) and the same model with 
no break in the CR source spectrum (red dashed line), are compared with available experimental data. In both cases we 
use DRAGON to model CR propagation and interactions (though almost identical results can be found with GALPROP). 
Here we use <I) = 450 MV to modulate both the Kolmogorov model and our Kraichnan reference models. <D = 300 MV 
was used only to match B/C ACE data which were taken in a very low activity solar phase. 



9 GV in the CR nuclei source spectra. The Kolmogorov model considered in adopts such a 
break. It is evident from Fig. |5(c)| as this is needed in order to reproduce the low energy tail of 
the observed proton spectrum which otherwise could not be fit for any choice of the modulation 
potential. It is important to notice, that this problem arises in all models with strong reacceler- 
ation x'A > 20 km ^ '. On the other hand our Kraichnan reference model requires a "modified" 
behavior of the diffusion coefficient at low energy (77 - -0.4 rather than t] - 1) which, however, 
may be motivated by independent physical arguments as discussed in Sec|4] 

From Fig.s |4(a)| the reader can see that while both the Kraichnan and Kolmogorov models 
reproduce the B/C equally well, the former model provides a significantly better description of 
the N/O ratio measured by HEAO-3 I30j] and CREAM llljLlsee Fig. [4(b]] i. Furthermore, what 
mostly favors our Kraichnan reference model are BESS ['33'], CAPRICE Ii47il and especially 
the PAMELA measurements of the p/p IT2I] and antiproton absolute spectrum llssll . Indeed, 
the discrepancy between low energy antiproton data and the prediction of the "conventional 
GALPROP model", which was akeady noted in |i6J, becomes more compelling due to the new 
PAMELA data, as shown in Fig. |5(b)| 

The comparison of our results with those of semi-analytical models is more difficult for ob- 
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Figure 5: The pip (panel a), p (panel b) and proton (panel c) computed with DRAGON are reported here use the same 
models and line notation used in Fig.s 4. The solar modulation potential used here is (D = 550 MV. In the (a) and (b) 
figures PAMELA 2010 data are taken from the recent paper B3l (see the note at the end of this paper). 
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vious reasons. One of the difficulties lies in the simplified gas and source distribution adopted in 
those models (see Sec.|2]i. We verified, however, that such difl'erences only affect the constraints 
to Do/z, with almost no eff'ect on the determination of 6. We also need to take into account 
that semi-analytical models (see e.g. |j2|,|3t]) assume diffusive reacceleration to take place only in 
the thin Galactic disk (whose height is Zd), while in the numerical models, as the one presented 
here, it takes place in the entire diffusion halo. Therefore, in order to compare the values of 
the Alfven velocity in those papers with those reported in the above it is necessary to perform 
a proper rescaling. This is approximatively given by (see e.g. Eq. (18) in [4]) va - v^/^ ^Zdlzt, 
with v^'^ being the Alfven velocity in the semi-analytical models and Zt the half scale height of 
the Galactic disk. 

In spite of these differences, and that CREAM and PAMELA data were not included in those 
analyses for chronological reasons, it is comforting that for low values of the convective velocity 
Vf ^ the preferred value of 6 estimated in |2, 3] is in remarkably good agreement with that 
found in this work: S - 0.45. Interestingly, the rescaled value of va determined in l|2i] is va - 10, 
for Vc - 0, which is also in good agreement with our results. It is important to notice that, 
similarly to what we did in our analysis, no break in the source spectral index was assumed in 

m. 

We remind the reader that in the above we always assumed Vc = as higher values of that 
parameter are not required to interpret CR nuclei and antiproton data. Models with a finite 
Vc, which may also allow to fit low energy data though with a different combination of some 
parameters (77, the modulation potential O, or any other), will be considered elsewhere. We 
already tested, however, that taking vv in a reasonable range of values do not affect significantly 
our constraints on most relevant diffusion coefficient parameters, namely 6 and va- Indeed, we 
verified that for various choices of the convective velocity, and of its vertical gradient, it is always 
possible to rescale the diffusion coefficient normalization Dq/z, so that both the B/C and the 
antiproton-to-proton ratio remain almost unaffected above 5 GeV/n, i.e. the energy range we 
considered in our analysis. Our tests confirm what is claimed in jst): namely, that the contribution 
of convection to the B/C energy slope is negligible, especially in the intermediate and high energy 
regions. 

6. Conclusions 

We used recent data on CR light nuclei and antiprotons to determine the conditions of prop- 
agation of high energy CRs in the Galaxy, exploiting our numerical code, DRAGON. In the 
framework of a diffusion-reacceleration model, we performed a thorough analysis of the agree- 
ment of our predictions with experimental information, aimed at constraining, in a statistical 
sense, the most important model parameters: Dq/zi, 6 and va- The amount and quality of data 
is enough to allow us to perform our analysis in a wide energy range, from 1 to 10^ GeV/n, 
and also to check the evolution of our results varying the minimal energy at which data are con- 
sidered. This is essential to reduce the uncertainties related to possibly unknown low energy 
physics, including solar modulation, and to disentangle the effects of reacceleration from those 
of diffusion. 

One of the most important results of this analysis is that light nuclei (especially B/C) data and 
antiproton data can fit into a unique, coherent diffusion-reacceleration model of propagation, as it 
can be read off Fig. [T] Indeed, the confidence regions obtained for E > 5 GeV/n (where only the 
effects of diffusion and re-acceleration matters), light nuclei and antiproton CL regions nicely 
overlap to produce combined constraints on Do/zt, S. While this was also shown in previous 
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works (which however did not exploit the new CREAM data), a combined statistical analysis 
of nuclear and antiproton data has been performed here for the first time. We showed that such 
an analysis allows to narrow significantly the allowed values of 6 and Dq/zi'- our constraints 
0.3 < S < 0.6 and 0.6 < Dq/zi < 1, as obtained at 95% C.L., are significantly more stringent that 
those previously determined in the related literature. Furthermore we found, for the first time, 
that only a relatively narrow range (10-20 km/s) of the Alfven velocity values are allowed. 
Even well below 5 GeV/n, we showed that it is possible to find effective models which, still 
fulfilling those constraints, allow to nicely reproduce all relevant data. 

We also found that the preferred values of the N/O and C/O ratios at injection are ~ 6 % and 
~ 75 % respectively. These results, and in particular the analysis of data with E,^in = 5 GeV/n, 
clearly favor a Kraichnan like CR diffusion (<5 = 0.5) respect to Kolmogorov (6 - 0.33). It is 
worth noticing that a relatively large value of 6, as that preferred by our analysis, would give rise 
to a too large CR anisotropy if our results are extrapolated to Ek » 10^ GeV/n (see e.g. |48] and 
ref.s therein). Our results, therefore, may call for some changes in the standard CR propagation 
scenario. 

While the effects of systematic uncertainties on fragmentation cross-sections are not studied 
in details in this work, recent results suggests that they should be smaller than the quoted 2cr 
statistical uncertainties on the transport parameters. Indeed we performed DRAGON runs using 
a different set of nuclear cross-sections as determined in ll40ll . finding the same B/C and N/O 
ratio found in the above within few %'s. 

Given that anyway nuclei data alone are able to provide constraints on Do/zr and 6, we use 
this information to establish a range for the maximal and minimal flux of antiprotons expected 
from CR interactions in the gas and still compatible with light nuclei observations within 95% 
CL. This range information can be used as a CR background in analyses aimed at constraining 
or finding some exotic signal in antiproton data. 

Forthcoming data from several running or scheduled experiments, as PAMELA (both for an- 
tiprotons and light nuclei), CREAM-II TRACER and AMS-02 HI which will 
measure both CR nuclei and p fluxes from hundreds MeV/n up to TeV/n, will soon allow tighter 
constraints. Especially AMS-02 is expected to provide very accurate data and, what is most 
relevant here, it will allow simultaneous and consistently calibrated measurements of several nu- 
clear species and antiprotons (as well as electrons and positrons which will also provide valuable 
complementary inputs). The AMS-02 potential to pinpoint the CR propagation was also recently 
showed in ll49[l where, however, the power of a combined analysis of CR nuclei and antiproton 
data was not discussed. 



Note added 

When this paper was in its final refereeing process, PAMELA collaboration published up- 
dated data of the ratio and p absolute spectrum |45]. The p/ p data differ very little from those 
we used in our statistical analysis (discussed in Sec. 3) which were taken from |12] (PAMELA 
2009). Therefore, their update should not affect significantly our constraints on the propagation 
parameters. In Fig. |5(a)| and |5(b)| we display the new PAMELA antiproton data. It is evident 
as our best fit model, which was determined using p PAMELA 2009 data, nicely matches also 
PAMELA 2010. 
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